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Abstract 



Quasielectrons and quasiholes in the fractional quantum Hall liquids 
obey fractional (including nontrivial mutual) exclusion statistics. Their 
statistics matrix can be determined from several possible state-counting 
scheme, involving different assumptions on statistical correlations. Ther- 
mal activation of quasiparticle pairs and thermodynamic properties of the 
fractional quantum Hall liquids near fillings \jm (m odd) at low tempera- 
ture are studied in the approximation of generalized ideal gas. The existence 
of hierarchical states in the fractional quantum Hall effect is shown to be 
a manifestation of the exclusonic nature of the relevant quasiparticles. For 
magnetic properties, a paramagnetism-diamagnetism transition appears to 
be possible at finite temperature. 

PACS numbers: 73.40.Hm, 05.30.-d 
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I. INTRODUCTION 



It is well-known that quantum statistics of a particle (or elementary excitation) plays a 
fundamental role in determining statistical or thermodynamic properties of a quantum many- 
body system. Bose-Einstein and Fermi-Dirac statistics are two well established ones, which 
are central to many familiar or novel phenomena involving many particles. For example, 
superfluidity or superconductivity is essentially due to Bose-Einstein condensation; and the 
stability of macroscopic matter is known to depend crucially on the Fermi-Dirac statistics of 
electrons. Since early days of quantum mechanics, an outstanding problem has been to search 
for a generalization of or even an interpolation between these two statistics. Mathematically, 
of course, there exist many possibilities. But as physicists we are interested in what are 
physically relevant, in the sense that the new statistics must be realized in physical systems 
existing in nature (or at least, in models which describe some interesting aspects of real 
physics) . 

In the last two decades or so, the interest in this search has become stronger and stronger 
in the study of lower dimensional condensed matter systems. For quite a while, it has been 
recognized that situations which interpolate between bosons and fermions may appear in 
one- [|1]J^] and two-dimensional || many-body systems, though no discussion of relevant 
statistical distributions until recently. 

By now there have been (at least) two distinct ways to define fractional statistics: 

1) by examining the change in phase of a multi-particle wave function due to the exchange 
of two identical particles, 

2) by counting the number of independent multi-particle quantum states to formulate a 
generalization of Pauli exclusion principle. 

For the usual quantum statistics, i.e. for bosons and fermions, the above two methods are 
equivalent to each other, in spite of their conceptual difference. For fractional statistics, 
however, they are generally inequivalent and we distinguish between the two definitions 
by calling them as "exchange statistics" and "exclusion statistics" respectively, and the 
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corresponding particles as "anyons" and " exclusons" . Fractional exchange statistics or 
anyons || have been first explored in the study of quasiparticles in two (space) dimensional 
systems, such as fractional quantum Hall (FQH) liquids and anyon superconductivity. 

Recently Haldane Q has formulated, by counting many-body states, a generalized Pauli 
exclusion principle in arbitrary spatial dimensions. Based on this idea, one of us (YSW) ]5[] 
have defined generalized ideal gas for particles obeying such fractional (including mutual) 
exclusion statistics, and have formulated its quantum statistical mechanics and thermody- 
namics. These new definitions are not merely mathematical construction; they have been 
shown to be realized as the exotic statistics obeyed by elementary excitations in certain one, 
two or higher dimensional strongly correlated systems [|J- [pjjj . In this paper, we discuss an 



important case: quasiparticles in the fractional quantum Hall (FQH) effect. 

By now it is well-known that the ground state of the two-dimensional electron gas in a 
strong perpendicular magnetic field with electron filling factor v = 1/m (m odd integer) is an 
incompressible quantum liquid |12[], and that it has two species of quasiparticle excitations, 



both of which are fractionally charged [0 (e?j_ = =)=l/m) anyons [|I^,|rj]]. However, the anyon 
approach is not very suitable for calculating low-temperature thermodynamic properties of 
the FQH liquids, since thermal activation of quasiparticle pairs is directly governed by the 
counting law for many-body states rather than the law for exchange phases. This is where 
the concept of exclusion statistics comes into play. The FQH quasiparticles are known to be 
strongly correlated. The key issue here is to clarify how strong correlation between quasi- 
particles manifests itself in the state-counting. Is possible not only fractional exclusion of 
single-particle states for identical quasiparticles, but also mutual exclusion of states between 
quasihole and quasielectron, which cannot be dealt with in the anyon approach. 

The many-body state counting for FQH quasiparticles is a subtle problem, whose study 
started with Haldane's paper in 1983 [|IIJ. Now we still do not have a full answer yet, ex- 
cept for the low-lying excited states of the FQH liquids, which however should be enough 
to account for low-temperature thermodynamic properties. There are several possible as- 
signments for the statistics matrix of FQH quasiparticles, involving different assumptions 



concerning the nature of the quasiparticles and their statistical correlations. (See below, Sec. 
Ill, for details.) Fortunately, these assignments of statistics matrix can be tested by numeri- 
cal simulation on small systems. In addition, one would also like to put these assignments to 
experimental tests, which needs comparing theoretical predictions with experimental data. 
As a first step towards this, we have calculated thermodynamic properties of the FQH liquids 
at low temperatures, based on the dilute gas approximation for thermally activated quasi- 
particles: When the quasiparticles are dilute, we may ignore Coulomb interactions between 
them, and apply the statistical thermodynamic formalism for generalized ideal gas given in 
ref. ||, which incorporates mutual statistics between different species of quasiparticles. It 
is hoped that this approximation could be improved in the future by including the effects of 
Coulomb interactions between quasiparticles. Our thermodynamic calculation is done with 
three different assignments of statistics matrix, with the hope that one day the experimental 
measurements of thermodynamic properties of FQH liquids might distinguish between them, 
providing information about the statistical correlations between FQH quasiparticles. 

This paper is organized as follows. We first review in Sec. II the state counting definitions 
of exclusion statistics, including mutual statistics between non-identical (quasi)particles. 
Then in Sec. Ill it is shown that quasielectrons and quasiholes in the FQH liquids with 
fillings 1/m (m odd) obey such fractional exclusion statistics, and the determination of their 
statistics matrix from various physical arguments or working hypotheses is reviewed, with 
several somewhat different outcomes. Using the statistical distribution |5j for generalized 
ideal gas, we study in Sec. IV and V, with analytic and numerical methods, thermal acti- 
vation of quasiparticle pairs near filling 1/m (m odd). Then we show in Sec. VI that the 
occurrence of new (hierarchical) incompressible states, corresponding to divergent pressure, 
at appropriate fillings at T = is a manifestation of fractional exclusion statistics, while 
at finite temperature the pressure of the system can never become divergent. In Sec. VII 
we compute low-temperature thermodynamic properties of FQH liquids, including magnetic 
properties. The final section VIII is devoted to conclusions and discussions. 
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II. EXCLUSION STATISTICS 



The definition of fractional exclusion statistics is directly based on state-counting, a basic 
concept in quantum statistical mechanics. It is well-known that bosons and fermions have 
different counting for many-body states, or different statistical weight W: The number of 
quantum states of N identical particles occupying a group of G states, for bosons or fermions 
respectively, is given by 

A simple generalization and interpolation is 

W = [G + {N-l){l- a )]\ 
N\ [G-aN- (1-a)]! 

with a = corresponding to bosons and a = 1 fermions. The physical meaning of this 
equation is the following: By assumption, the statistical weight remains to be a single 
combinatoric number, so one can count the states by thinking of the particles effectively as 
bosons, with the effective number of available single-particle states being linearly dependent 
on the particle number. 

G% = G-a{N-l). (2.3) 

Obviously, for genuine bosons , Gfh is independent of the particle number. In all other 
cases, G^jf is linearly dependent on the particle number. This is the defining feature of 
the fractional exclusion statistics. The statistics parameter a tells us, on the average, how 
many single-particle states that a particle can exclude others to occupy. For a ^ 1, this 
generalizes the Pauli exclusion principle for one species. 

It is easy to generalize this state counting to more than one species: 

w = Tf [fr + jVi-l-S^W-frj)]! (2A] 

n * (iVj)! c, i v ; n, ; (.Y,-rt; ; );! • 1 • ] 

Here Gi is the number of states when the system consists of only a single particle. By 
definition, the diagonal otu is the "self-exclusion" statistics of species i, while the non- 
diagonal cx.ij (for i j) is the mutual-exclusion statistics. Note that a^, which Haldane 0] 



called statistical interactions, may be asymmetric in i and j. The interpretation is similar 
to the one species case. The number of available single-particle states for species i, in the 
presence of other particles, is again linearly dependent on particle numbers of all species: 

G% = G ^-J2 a dN j -S ji ). (2.5) 

3 

We note that as a generalized Pauli exclusion principle, eq. (|2.2j ) or (|2.4j ) implies strong 



correlations between the particles, and does not give rise to the same state counting as the 
old generalization suggested in ref . fl6| , in which particles independently fill a fixed number 



of single-particle states with the constraint that at most n particles are allowed in one and 
the same state. 

Some remarks on exclusion statistics are in order: 

1) This definition of exclusion statistics is independent of spatial dimensionality of the 
system, in contrast to the exchange statistics of anyons which has a connection to the braid 
group, making sense only in two spatial dimensions ||17|| . 



2) In contrast to anyons, there is no periodicity in exclusion statistics parameter a, so it 
makes sense to consider the cases with a > 1 or even a > 2. 

3) The state-counting definition of exclusion statistics naturally allows mutual statistics 
from the beginning, implying that exclusion may occur between states of different species, 
a completely new situation we have not been faced before in statistical mechanics. 



III. EXCLUSION STATISTICS FOR FQH QUASIPARTICLES 

There are two kinds of quasiparticles in the Laughlin 1/m-liquid: quasiholes labeled by 
— and quasiparticles labeled by +. In this paper we treat them as two distinct species 
and demonstrate that their many-body states obey the counting law given by ( |2.5| ), with 
an appropriate 2x2 statistics matrix aij with i, j = +, — . The statistics matrix depends 
on the nature of FQH quasiparticles and their correlations. Several scenarios are possible 
in this regard. In this section, we are going to discuss four possible scenarios for FQH 
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quasiparticles that have appeared in the literature: 1) anyons in the lowest Landau level, 
2) bosonic vortices, 3) composite fermions, 4) correlated vortices or composite fermions. 
They differ in the assumption of whether certain correlations, such as hard-core constraints, 
exist between the quasiparticles or not, leading to subtle difference in statistics matrix. It is 
remarkable that the statistics matrix can be subject to numerical test for small systems on 
a sphere. We are not going to talk about the details, but will briefly summarize the status 
of such numerical tests and quote relevant references when appropriate. 



A. Anyons in the lowest Landau level 

In determining exclusion statistics of FQH quasiparticles, let us first try to explore the 
fact that they are fractionally charged anyons. (Though later we will see that the picture of 
non-interacting anyons is not very suitable for calculating thermal activation of quasiparticle 
pairs.) 

Good trial electron wave functions for states with quasiparticles in the 1/m FQH liq- 
uid were first proposed by Laughlin ]TJ|. In these wave functions the coordinates of the 
quasiparticles appear as parameters (or collective coordinates). If one moves very slowly the 
coordinates of one quasiparticle, say a quasihole, around a closed loop in the FQH liquid, 
the electron wave function acquires a Berry phase, which can be interpreted as the phase 
due to the motion of the quasiparticle traveling along the loop. As shown in ref. ||14j| , the 
Berry phase is always proportional to the number of electrons enclosed in the loop. If the 
loop encloses none of other quasiparticles, the Berry phase is the same as that for a charge 
in a magnetic field, with electrons acting as quantized sources of "flux". Thus, a quasi- 
particle sees the electrons just like an electron sees the external magnetic field. When the 
loop encloses another quasiparticle, say a quasihole, the change in the Berry phase is due 
to a deficit in the number of enclosed electrons caused by the enclosed quasihole, and it is 
attributed to the exchange phase of the two quasiholes, showing that they are anyons with 
fractional exchange statistics |H| . Combining the two results, one is led to a simple picture 
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that the FQH quasiparticles are anyons in the lowest Landau level of a fictitious magnetic 
field, whose strength is determined by the density of electrons. Indeed, the wave function 
for the FQH quasiparticles suggested by Halperin |18[ are such as if the quasiparticles are 
in the lowest Landau level. 

Now let us count the states of N non-interacting anyons (of one species) in the lowest 
Landau level. Though not for all levels and all states, a number of exact solutions for anyons 



in a magnetic field have been known JL{|. Among them, fortunately, are the complete set 
of solutions for all anyons in the lowest Landau level (if the number of anyons is less than 
the Landau degeneracy). The total (ground state) energy turns out to be the sum of the 
cyclotron energy of individual particles, independent of the exchange statistics of anyons: 

E = Ne c /2 . (3.1) 

To count the states, we consider anyons in a circular disk with a fixed size. In the 
symmetric gauge, besides the usual Gaussian factor, the many-anyon wave function is known 
to be of the form (with Z{ the complex coordinates of electrons): 

* = J[(z i -z j ) ' v -$(z l ,...,z N ), (3.2) 

i<j 

with anyon statistics < 9 < 2n. But now in the lowest Landau level, the function $ is a 
symmetric polynomial of (z±, . . . , z^). The state counting can be easily done by looking at 
the symmetric polynomial $ (counting as bosons). However, the fixed-size condition requires 
a fixed highest angular momentum, or a fixed highest power of a single variable Zj in the 
wave function On the other hand, the Jastrow-type prefactor nj<j(zj — Zj) e ^ n , implying 
non- vanishing relative angular momenta between anyons, takes away some powers of Z{ and 
reduces the degree of the polynomial $. Alternatively, an anyon can see the statistical flux of 
other anyons, which in the present case is opposite to the external magnetic flux. Therefore, 
in the boson counting, with size and external flux fixed, the effective Landau degeneracy is 
determined by the external magnetic flux less the anyon statistical flux {9/n){N — 1): 



at a 

G% = £l±-^N-l). (3.3) 
IJ m 7r 

Hence eq. (|2.2j) applies, with the single anyon Landau degeneracy G = N^/m, and the 
exclusion statistics for anyons in the lowest Landau level can be read off from eq. (|2.3|) : 



a = 0/tt. (3.4) 
In Ref. [[H], the exchange statistics has been shown to be 6L = ix/m for quasiholes, and 



9 + = —7i/m for quasielectrons. Thus, eq. ( |3.4j ) leads to the following diagonal exclusion 
statistics for quasiparticles: 

1 (—1/m (soft — core), 

a = — , a ++ = < (3-5) 

^ 2 — 1/m (hard — core). 



m 



For the case of quasielectron, we note that there are two possibilities: Since conceptually 
exchange statistics is an angular parameter, defined only up to a period of 2n, 9 = —n/m is 
equivalent to 9 = (2 — l/m)ir. However, exclusion statistics is always unambiguously defined 
and non-periodic at all; i.e. the exclusion effects with a = —1/m and with a = 2 — 1/m 



are very different. So when one wants to apply eq. ( |3.4j ), he or she has to choose between 
the possible two values of 9. One may notice that the wave function ( |3.2| ) with 9 = —n/m 
is singular at z$ = Zj. Based on the braid group, one may argue that the many-anyon 
wave function Q2.1|) should vanish as two anyons approach each other, thus preferring 9/tt = 
2 — 1/m over —1/m in the Jastrow-type prefactor and therefore a = 2 — 1/m over —1/m for 
the exclusion statistics for quasielectrons. Obviously, the former value of a leads to stronger 
exclusion between quasielectrons, as if they have "hard-core". 

Whether the quasielectrons really satisfy the "hard-core constraint" can be tested by 
numerical experiments. Such numerical experiments have been done by three groups, p0| . 



21 1 and |p2 |, for electrons with Coulomb interactions on a sphere in the field of a monopole 
at its center. Their results unambiguously support the exclusion statistics ( |3.5| ) with the 
"hard-core" value for quasielectrons. 
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B. Bosonic Vortex Scheme 



The above scenario for FQH quasiparticles as anyons in the lowest landau level has the 
disadvantage that it tells us nothing about the mutual (or non-diagonal) statistics between 
quasihole and quasielectron, which is important for studying thermal activation of quasi- 
particle pairs. So we need other, more direct ways to count states for quasiparticles in the 
Laughlin 1/m-liquid. 

A fundamental relation which will play a key role in state-counting is the "total-flux" 
constraint 

= eBV/hc = mN e + AL - N+ , (3.6) 

between the electron numbers N e and quasiparticle numbers iV_ and N + . The basic idea 
behind this relation is the following observation |L2| : To generate a quasiparticle in the 
incompressible Laughlin liquid, one may pierce the droplet by an infinitely thin solenoid and 
slowly turn on magnetic flux inside it. When the flux reaches a flux quantum, a quasiparticle 
will be formed around the solenoid; whether it is a quasihole or quasielectron depends on 
the direction of the solenoid flux (parallel or anti-parallel to the external magnetic field). 

Further state counting relies on the assumptions on statistical correlations of quasiparti- 
cles. There are different counting schemes based on the bosonic vortex picture, the composite 
fermion picture and variation of both. Let us consider them in turn. 

The bosonic vortex scheme is based on the picture that the FQHE quasiparticles are 
vortex-like excitations in the incompressible planar quantum liquid, and for a fixed number 
of excitations we count their states as if they are bosons. Assuming only the minimal 
(quantized) circulation, there are two possible orientations for vortex circulation on the 
plane, corresponding to quasihole and quasielectron respectively. So what is essential to this 
counting scheme is to determine the number of available states for each species of vortices. 



This can be inferred from an observation by Haldane and Wu [23] that for vortices in a 



planar quantum liquid, their core X- and F-coordinates do not commute with each other, 
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as if they were the guiding-center coordinates for a charged particle in a magnetic field, 
with fluid particles (i.e. electrons in the present case) as sources of quantized flux. Thus, 
the number of available states for vortex-like excitations is essentially determined by the 
"Landau degeneracy" of this fictitious magnetic field, or the number of electrons: 

G% t _ = N e , G?} fi+ = N e . (3.7) 

To derive the exclusion statistics for FQH quasiparticles, one needs to fix the external 
flux Nfj,. So let us express G e ff tT in terms of N<f, by eliminating N e from these equations 



with the help of the constraint (|3.6| ). Then we obtain 



Gf} f _ = -N^--N_ + -N, 
JJ ' m m m 



G% + = -Nf - -AL - (--)JV+ . (3.8) 

The first term on the right side gives the single-quasiparticle degeneracy in terms of the 
external flux: 

G + = G_ = (l/m)A^, (3.9) 

so the proportionality constant 1/m is identified as the fractional charge (absolute value) of 
the quasiparticles. And the coefficients of N T give the statistics matrix: 

a ++ = —1/m, ct+_ = 1/m, 

(3.10) 

«_ + = —1/m, a = 1/m . 

This result was first derived by Haldane 01]. 

Comparing with eq. ( |3.5| ), we note that a ++ = —1/m here corresponds to soft-core 
quasielectrons. For hard-core quasielectrons, eq. ( p.7| ) should be replaced by 

G$j ft _ = ^e, = N e - 2(N + - 1), (3.11) 

with the second term in Gfh + representing the exclusion effects due to the hard core of 
quasielectrons. In the presence of an external magnetic field, the two orientations of vortex 



1 „ 1 „ , 1 
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circulation are not equivalent, so there is an asymmetry between quasiholes and quasielec- 
trons. Physically, the hard-core nature of quasielectrons may be due to electron number 
surplus in the core of quasielectrons. This leads to [^4]] 

G^l_ = -N^--N. + -N + , 
JJ ' m m m 

= -JV, - -JV_ - (2 - -)N + , (3.12) 
jj,t m m m 

resulting in a statistics matrix somewhat different from eq. ( ft. ICQ : 

a++ = 2 — 1/m, a+- = 1/m, 

(3.13) 

a |_ = —1/m, a = 1/m . 

We note that now the diagonal statistics a ++ in eq. ( ft,13|) agree with that of hard-core 
quasielectrons in eq. ( ft.5| ). However, here it has been possible to demonstrate nontrivial 
mutual statistics between quasihole and quasielectron. This means that with fixed, the 
presence of quasielectrons will affect the number of "available" states for quasiholes and vice 
versa. 



C. Composite Fermion Scheme 

The central idea of the composite fermion approach is that the FQH state of electrons 
in a physical magnetic field can be explained as the IQH state of composite fermions in an 



effective magnetic field [p5 1 . Imagine an adiabatic process in which we somehow collect 2p (p 
an integer) flux quanta to each electron to form an electron-flux composite. The additional 
Aharonov-Bohm phase, due to the attached flux, associated with moving one composite 
around another is e l2p7T = 1. So the statistics of the composite remains to be the same as 
the electron, motivating the name of composite fermion. These composites are now moving 
in a reduced magnetic field B e ff = B — 2n(2pp), where p is the density of electrons, which 
is the same as the density of composite fermions. (Recall that in our convention, the unit 
of flux is 2tt.) The filling factor for composite fermions then increases to v e ffi given by 
v~ff = (B — 4:7rpp)/27ip = v~ x — 2p. For u e ff = n (n an integer), we have 
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Thus, fractional Hall systems with v = n/(2pn + 1) may be adiabatically changed into an 
integer Hall system with filling factor n, as was also emphasized by Greiter and Wilczek 
26fl . Note that this argument gives us more than we had hoped for. The case we wanted to 



understand, with z/ _1 = odd, is obtained for n = 1. 

Let us do state counting for the state with n — 1, or v — l/(2p + 1). The Landau 
degeneracy for the composite fermion in the residue magnetic field is given by the effective 
flux 

N 4>ieff = N <f> -2pN e , (3.15) 

while the number of excitations are determined by 

N Mf = N e + N--N+, (3.16) 

Eliminating N^eff from these two equations, we recover the same constraint ( |3.6|) as before 
with m = 2p + 1 . 

The number of available single-particle states for unit-charged composite-fermion exci- 
tations is obviously 

G eff , T = N^ ff -(N T -l). (3.17) 

Here we have assumed that the magnetic field is so strong that we can ignore the possibility 
for quasielectron to fill Landau levels higher than the lowest available one. 

To derive the true charge and statistics of the quasiparticle excitations, we need to express 
N^^ff in eq. ( |3.17| ) in terms of the external N^, resulting in |fZ|[ 



G@ /f _ = -N* - -AL - (1 - -)iV 4 
■ ■ m m m 



1 „ 1 „ „ 1 



G% + = -N+ + (1 - -)AL - (2 - -)N + . (3.18) 

Here we have used eqs. ( |3.15| ) and ( |3.16| ). The coefficient of recovers the fractional 
charge (absolute value) 1/m for the quasiparticles; see eq. ( J3.9|) . From the coefficient of the 
other terms one reads off the exclusion statistics: 
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1 „ L„ 1 



a ++ = 2 — 1/m, a+- = —1 + 1/m, 

(3.19) 

a_ + = 1 — 1/m, a = 1/m . 



D. Correlated Vortex or Projected Composite Fermion Scheme 



Comparing eq. (|3.19[) with eq. (|3.13|) , we see that the composite fermion scheme leads 



to the same diagonal statistics both for quasiholes and for quasielectrons as the bosonic 
vortex scheme, with quasielectrons being automatically hard-core. Thus it is not surprising 
that the two schemes give the same prediction about the occurrence of hierarchical states at 
T = 0, since only one species of quasiparticles exist at T = when the filling factor deviates 
from the magic 1/m, so that only diagonal statistics is relevant. 

However, mutual statistics in eq. ( |3.19| ) and eq. ( |3.13| ) obtained from the above two 



schemes are obviously different. Which is correct? Or neither is correct? To decide, one 
needs to study situations in which both species of quasiparticles coexist at the same time. 
This problem has been numerically studied in ref. |27] (see also |29|). It turns out that 



neither of the mutual statistics given in eq. (|3.13|) and eq. ( 3.19D is correct. The correct 
statistics matrix, for low-lying excitations, turns out to be 

a ++ = 2 — 1/m, a i = —2 + 1/m, 

(3.20) 

ct_ + = 2 — 1/m, a = 1/m. 

In the bosonic vortex scheme, this can be obtained by incorporating certain amount of 
mutual exclusion (or inclusion) between vortices and anti-vortices in eq. ( 3.1 1|) as follows: 



Geff,- = Ne~ 2iV 4 



G% t+ = N e + 2M - 2(N + - 1) . (3.21) 



We call this scheme as the correlated vortex scheme, since in ref. |2^J it has been shown 
that this modification is due to the necessity of inserting some "hard-core" Jastrow factor 
between quasihole and quasielectron in the quasiparticle wave functions, which represents 
correlations of a new type between the vortex and the anti-vortex. 
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To reproduce the statistics ( p.20|) in the composite fermion scheme, one needs to modify 
eq. ( ETT7D to |28],|29| 



G? e//l _ = N Mf - (JV_ - 1) - JV + , 
G effl+ = N Mf + N- - {N+ - 1) . 



(3.22) 



The mutual exclusion added in these equations between composite fermionic holes and com- 
posite fermionic electrons can be interpreted as a consequence of the projected composite 
fermion scheme: In the composite fermion transformation quasielectron states involve wave 
functions in the second Landau level, so it is necessary to project the wave functions down 
to the lowest Landau level to obtain the correct many-electron wave functions. (For details, 
see ref. [SOfl.) Indeed, the state counting resulting from the above equation has been checked 
28|,29] to be indeed in agreement with the numerical data given in ref. [pCj, which verifies 



the necessity for the projected composite fermion scheme. 

We note that in either scheme, the mutual statistics between quasihole and quasiparticle 
are anti-symmetric rather than symmetric. 

In summary, the bosonic vortex scheme ( 3.1 lj ) and the unprojected composite fermion 
scheme ( |3.17| ) lead to different mutual statistics ( |3.13| ) and ( ft.l9| ) respectively. But the 
correlated vortex scheme ( 3.21]) and the projected composite fermion scheme ( |3.22|) lead to 
the same statistics matrix (|3.20| ). Numerical data favor the latter. But we feel that there is 
no harm to leave these possibilities open to experimental tests. In the following, we calculate 
thermodynamic properties of FQH liquids with the three statistics matrices, with the hope 
that someday experiments might be able to distinguish between them. 



IV. THERMAL ACTIVATION OF FQH QUASIPARTICLE PAIRS 

It is well-known that the low-temperature thermodynamics of a many-body system is 
determined by its low-lying excited states above the ground state. For FQH liquids at 
hand, our fundamental assumption is that their low-lying excited states are dominated by 
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weakly coupled quasiparticles. Thus their low-temperature thermodynamic properties are 
determined by thermal activation of FQH quasiparticle pairs. At low temperatures, when 
the activated pairs are not very dense, one may ignore their interaction energies. Then the 
densities p± of the excitations should be determined by the laws for generalized ideal gas 
with two species, in which the following two conditions are satisfied: 

1) The state-counting (|2.4| ) for statistical weight W is applicable. 

2) The total energy (eigenvalue) is always of the form of a simple sum, in which the i-th 
term is linear in the particle number Ni m . 

E = Y J N l e i , (4.1) 

i 

with Si identified as the energy of a quasihole {i = — ) or a quasielectron {i = +). Though 
this condition ( ^4. 1|) is very natural for weakly interacting FQH quasiparticles, we note it 
is not compatible with non-interacting anyons, except for only one species of anyons all 
in the lowest Landau level (see eq. ( |3.1| )). This problem does not exist in the theoretical 
framework of exclusion statistics: The condition ( |4.1| ) is compatible with free exclusons, 
as exemplified 0- [|TTJ in one-dimensional exactly solvable many-body models such as the 



^-function repulsive boson gas [Q] and the Calegero- Sutherland model 0. This is one of 
the main theoretical advantages of exclusion statistics over exchange (or anyon) statistics 
in dealing with statistical mechanics. (Moreover, the anyon picture can not deal with more 
than one species, so it is not suitable for studying thermal activation, which involves both 
quasielectrons and quasiholes and is expected to be a good place to look for the effects of 
mutual statistics, with increasing density of activated pairs.) 

With theses assumptions, now we are able to derive quantum statistical mechanics of 
FQH quasiparticles. Consider a grand canonical ensemble at temperature T and with chem- 
ical potential /i« for species i = +, — , whose partition function is given by (with k the 
Boltzmann constant) 

Z = ]T WdN,}) exp{ £ - e t )/kT} . (4.2) 

W;} i = +- 
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As usual, we expect that for very large Ni, the summand has a very sharp peak around 
the set of most-probable (or mean) particle numbers {Ni}. Using the Stirling formula and 
introducing the average "occupation number per state" defined by rii = Ni/Gi, from the 
maximum condition 

A [log W+ G jnj ( H - e 3 )/kT] = , (4.3) 
1 i=+- 

one obtains the equations determining the most-probable distribution of rii 

E (Vi(r) + ^ = l. ( 4 - 4 ) 
j=+- 

with Wi(T) being determined by the functional equations 

w + a++ (l +w + ) 1 - a++ ( W ~ Y~ + = e ^+-^)/kT 

w_---(l + w^ 1 - a --(^—y + ' = e< e — »~y kT . (4.5) 

From eqs. (fOJ) and ( [4.5[ ), n± are expressed in terms of w± as 

p+(T) W- + a — a+_ 



n + {T) 



po (w+ + a ++ )(w- + a — ) - a+_a_+ 



Po (w + + a ++ ){w^ + a ) — a+_a__ 



-+ 



where po = G±/V, and p+(T) and P-(T) are the density of quasielectrons and quasiholes 
respectively. The ratio -R(X') of the numbers of quasielectrons and quasiholes is given by 

R(T) = ^ = W -l a — a +- . (4.7) 
n_(i J w+ + a++ — a_+ 

According to charge conservation, only quasielectron-quasihole pairs are thermally ac- 
tivated, since they have opposite charges. Thus, iV + — AL is independent of temperature. 
Then the total-flux constraint ( |3.6| ) implies that 



n+ (T) -n_(T) = m6, (4.8) 

where 5 = m[y — uq), (uq = 1/m). Thus 8/m gives the deviation of the filling from 1/m. 
Then from eq. ( |4.6|) one has 
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m6 = + + . (4.9) 

(w + + a ++ )(u>_ + a ) — a + _a_ + 

Charge conservation also requires that the chemical potentials for the two species should 
satisfy p + + /i_ = 0. Multiplying the two equations in ( |4.5| ), and using the above constraints, 
one can derive a polynomial equation 

w+ <*+++<*+- w _<x-++<*—(i + w+ y- a ++- a +- (i + w _) l - a -+~ a -- = e ( £ ++ e -)/ fcT . (4.io) 

Once w + and w;„ are determined from eqs. ( |4.9| ) and ( |4 . 1 U| ) , the T-dependent densities 



of both species p+ = n + p and p_ = ri-p are given by eq. (j4~q ). The mutual-statistics- 



dependent thermodynamic potential Q = — kT log Z and entropy S are 

Q = _pv = _& T £ GilQg Po + Pi-^=+,-^Pj . (4 n) 

i=+,- Po ~ T,j=+,~ a ijPj 

f = E ^{n^ + log P ° + P V E ^- a ^ }. (4.12) 
Further the total entropy is written as 5* = J2% ^i s i with 

T = [1 + - E «of ] log[l + - 

ft Pi j Pi j Po 

-log^-(^-E^-) log(l-E^-)- (4-13) 

PO Pi j Pi j Po 

Other thermodynamic functions, such as specific heat and magnetization per unit area, 
follow straightforwardly. For example, magnetization per unit area is given by 

Here fj,± = de±/dB; and we have assumed the independence of fi±. Hopefully, when T is of 
order of e± or higher, the ay-dependent second term may give an appreciable contribution. 
Note that the thermodynamic properties at the two sides of electron filling uq = N e /N^ = 
1/m are not symmetric, due to asymmetry of quasielectron and quasihole in self-exclusion 
and mutual-exclusion statistics. 
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V. EXPLICIT SOLUTIONS NEAR u = l/M 



We present some explicit formulas for thermal activation of FQH quasiparticle pairs for 
three statistics matrices ( 3.13| ), ( 3.19 ) and (|3.2U| ), which were derived from three different 



counting schemes in Sec. III. 

A. Bosonic Vortex Picture 

From the statistics matrix (|3.13|) , eq. (|4.6| ) determines the occupation number of quasi- 
electrons and of quasiholes to be 

n = (5 x) 

+ (w+ + 2 - l/m)0_ + l/m) + l/m 2 ' V ' ; 



and 



(w + + 2 - 1 /m) (w_ + l/m) + l/m 2 



where w + and u>_ can be obtained by solving eqs. ( |4.9| ) and ( [4.1(J| ), which are explicitly 



m5w + w„ + (5 + l)w+ + [(2m - 1)5 - l]w- + 2 + 25 = 0; (5.3) 

wl(l + w + )-\l + «,_) = e A / fcT = /(T), (5.4) 

where A = e + + is the pair excitation gap. 

The last two equations can not be solved analytically, but numerical solution is possible. 
See Figures |l] (a) and || (a) for a two-dimensional plot of n + (T,5) and n_(T, 8), obtained 
numerically for near l/m = 1/3. 

B. Composite fermion scheme 

In this scheme, the statistics matrix is given by eq. (|3.19|). Then eq. (|4.(j|) is explicitly 



n = w - + 1 (55) 

+ (w+ + 2 - l/m)(w- + l/m) + (1 - l/m) 2 ' v ' ; 
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and 

n ~ (w+ + 2 - l/m)(u;_ + 1/m) + (1 - 1/m) 2 ' 1 ' j 

with w + and w_ determined by 

«>+«;_ = e A/fcT s /(T), (5.7) 

mSw + w_ + (5 + l)w+ + [(2m - 1)5 - + m<5 = 0. (5.8) 
In this case, the solution in analytic form is available: We obtain explicitly 
w + (T) = yf^Y ) {- ml >\l( T ) + 1] ± V^W) + I] 2 - 4(5 + l)[(2m - 1)<S - 1]/(T)}, 
- ^ • (5-9) 



The upper sign is for 5 > and the lower sign for 5 < 0. At T = 0, eq. (|5.9| ) indeed yields 
n + = m5, n_ = 0, for 5 > 0, and similarly n + = 0,n_ = — m8, for 5 < 0, as expected from 
eq. O. 

Numerical results for a two-dimensional plot of n+(T, 5) and n_(T, 5) for filling factors 
near 1/m = 1/3 are shown in Figures [l| (b) and (b). 

C. Correlated Vortex or Projected Composite Fermion Picture 

For the statistics matrix ( |3.13| ), eq. ( |4.G| ) is explicitly 

„ = w - + 2 ( 5 i ) 

+ (w + + 2-l/m)(u;_ + l/m) + (2-l/m) 2 ' V ' ; 

and 

n _ = ^± (5 id 

{w + + 2- l/m)(w_ + 1/m) + (2 - 1/m) 2 ' V ' ' 

with w + and w_ satisfying 

m6w + w- + (5 + l)w+ + [(2m - 1)5 - l]ty_ + 2(2m - 1)5 - 2 = 0, (5.12) 
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w 2 _(l + w+)(l + = f(T). (5.13) 

Again, analytic solution is impossible, but a numerical two-dimensional plot for n + (T, S) 
and n_(T, 5) is shown, respectively, in Figures [1] (c) and |2] (c), for filling factors near 1/m = 
1/3. 



D. Low- Temperature Asymptotics 

As application of the above explicit formulas, let us discuss the low-temperature asymp- 
totics of the density of activated pairs. For simplicity, we consider the case with exactly 
v — 1/m, or 5 — 0. It is easy to check that in either of the schemes, we have n + (T) = n_(T). 
At very low temperatures, f(T) is very large, so we have w± ~ exp{A/2fcT}. This leads to 

p±(T)^p exp{A/2kT}, (5.14) 

with the prefactor p = G±/V = (l/m)^/]/, proportional to the (fractional) quasiparticle 
charge. 

This is in complete agreement with the standard Boltzmann behavior characteristic of 
thermal activation across a finite gap. Note that this behavior is independent of the statistics 
matrix. Thus to look for the effects of fractional exclusion statistics, the temperature should 
be higher than this asymptotic region. 



VI. EMERGENCE OF HIERARCHICAL STATES 



From the general equation of state ( |4.11|) , one can easily see that the pressure P will di- 



verges, when one of the denominators in the right side becomes zero, i.e. when the excitation 
densities satisfy 

aijUjiT) = 1 (i = + or -). (6.1) 

3=+,- 
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This corresponds to the situation in which one of the G^h^ vanishes (see eq. (|2.5|) ), so that 
there is no available quasihole or quasielectron states for additional thermally activated pair 
to occupy. Using eq. (|4.4|), the condition (|6.2|) is reduced to 



w + n + = 0, or = 0. (6.2) 

At zero temperature in the ground state with the filling factor v near the magic z/ = 1/m 
(m odd), there is no thermally activated pair, so one of the n± vanishes, depending on 
whether v is greater or smaller than 1/m. There are three cases that the condition (|6.2|) is 
satisfied: 

(i) n + = ri- = 0: This is the case with v = uq, which is just the original Laughlin's 1/m 
incompressible state. 

(ii) n_ = and n + > 0: According to the relation ( |4.8|) , this is the case with v > u . Then 
w + = 0, and eq. ( |4.10| ) requires w_ = oo. So from eq. (fOl), one obtains n + = l/a ++ . 



For a ++ = 2 — 1/m, it leads to v = 2/ (2m — 1). This gives rise to a new incompressible 
state in which quasiholes are absent and quasielectrons fill up all possible states. In the 
bosonic vortex scheme, this new quantum Hall state is called the first hierarchical state 
IRIS], in which the hard-core quasielectrons in the Laughlin 1/m-state are condensed to 



form a new incompressible liquid. If quasielectrons did not have hard core, they would not 
be able to form the new incompressible hierarchical state. In the composite fermion scheme, 
the original Laughlin 1/m-state corresponds to n = 1 in eq. (|3.14j ), and is interpreted as 



complete filling of the lowest Landau level in the residue magnetic field, while the new state 
corresponds to n = 2 with the second Landau level completely filled by composite fermions. 
(iii) n + = and n_ > 0: According to the relation ([4.8|) , this is the case with v < uq. Then 
w_ = 0, and eq. ( [1.10| ) requires w + = oo. So from eq. ([4.6|), one obtains = 1/a For 



a = 1/m, it leads to v — 1 — 1/m state, in conflict with the condition v < u . So unlike 

quasielectrons, the quasiholes in the Laughlin 1/m-liquid can not condense to form a new 
incompressible state. Numerical data presented in ref. |2D[ confirms this conclusion. 



We note that in the above discussions for T = 0, mutual statistics is irrelevant. 
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What will happen at finite T? Suppose that v > uq. Due to thermal activation of 
quasiparticle pairs, both n±(T) > 0, so there is an additional contribution in eq. ( p . 1| ) 
from thermally activated quasiholes, which is T-dependent. One may wonder if the effect 
of mutual statistics would lead to formation of new incompressible (hierarchical) states at 
T-dependent filling factors. At finite T, the condition ( p.'2\) is reduced to w + (T) = or 
W-(T) = 0. It can be shown, case by case for the three statistics matrices ( [3.13|) , (|3.19|) 
and (|3.20| ), that in either scheme, the pressure can not diverge at any filling factor. For 
example, in the composite fermion scheme, eq. ( |5.7|) with finite T implies that if w+(T) = 0, 
then w-(T) = oo, which in turn leads to n_(T) = 0, in accordance to eq. (|5.6| ). But this 
is impossible at finite T due to thermal activation. Similarly iw_(T) = would lead to 
n + (T) = 0, a contradiction too. The same is true in the other two schemes. 

The fact that the pressure P never diverges at finite T means that the quantum Hall 
transitions due to a divergent P are a quantum phase transition at zero temperature. 



Under the assumption that quasiparticles dominate the low-lying excitation spectrum 
of the FQH liquids, the knowledge of thermal activation of quasiparticle pairs allows us 
to calculate thermodynamic observables of the FQH liquids. In this section, we will show 
numerical results with three possible statistics matrices ( |3.13|) (|3.19|) and (|3.20| ). 

In Figure |3], we show the temperature dependence of the thermonynamic potential pV 
for several different filling factors. 

The average energy density h(T) of the quasielectron-quasihole pairs is given as 



VII. THERMODYNAMIC OBSERVABLES 



h(T) = e + [(p + (T) - P+ (0)] + e-[{P-{T) - p_(0)]. 



(7.1) 



The specific heat C V (T,5) is also evaluated numerically using the expression 



dh(T) 



-kp E 



2 dn - 



(7.2) 



V 



dT 



BE ' 
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where E = A/kT. In Figure ^, we have shown the curves of C V (T) at different values 
of S for the three statistics matrix. At very low temperatures, C v vanishes exponentially: 
C V {T) ~ exp{— A/kT}, due to a finite activation gap. At higher temperatures, C v increases 
due to quasiparticle pair creation. When the temperature approaches the order of magni- 
tude of the gap energy A, the difference between different mutual statistics becomes more 
and more apparent. The phenomenon that C V (T) decreases after reaching a maximum is a 
manifestation of exclusion statistics due to the saturation of the available states for quasi- 
particles. The position of the maximum of the C V (T) curve is statistics-matrix dependent. 
Whether our approximation is still valid at this temperature or not is a question we cannot 
answer in our approach. We leave it to experiments. However, we are sure that the van- 
ishing of C v at very high temperatures should not be trusted, since our approximation of 
restricting quasielectrons to the lowest available Landau level certainly breaks down. 

Now let us consider the magnetic response of the system. The magnetization per unit 
area due to quasiparticle excitations has been given in ref. as 



M 



\ dB ) 



T,V <fH /kT 



^(-,f p , + ^L la l±^y (7 .3) 



»=+,- 



?0 Wi 

Here <po = hc/e. In the derivation, we have effectively defined fj,± = de±/dB. In comparison 
with experimental data, the effective magnetic moment, of a single quasiparticle could 
be either treated as phenomenological parameters or derived from some microscopic model. 
In the following, we will assume \i± to be ^-independent. 

This equation can be viewed as a generalization of the case with one species of anyons 
in the lowest Landau level. [[H],[5] The first term, corresponding to the usual de Haas-van 
Alphen term, is of the same form in the two cases. However, there are important differences 
between the one-species anyon case and the present situation: For the former case, the 
density of anyons, p, is fixed either when B or T varies (no thermal activation of anyon 
pairs). The resulting susceptibility has a simple analytic expression and vanishes as T — ■> 0. 
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|| However, this is not true in the present case, where the densities of quasielectrons and 
quasiholes are both B- and T-dependent. Even at zero temperature, the magnetization has 
a independence through that of p±: 

M(T = 0)=p^ P± (B,T = 0)=p^p n ± (B,T = 0). (7.4) 

(The subscript ± depends on whether u — l/m is positive or negative.) This implies that the 
zero-temperature susceptibility is non-zero, in contrast to the case of anyons of one species: 

Our general expression for the finite-temperature susceptibility reads 

'dM\ 



X 



where Xo is a de Haas-van Alphin-like term and Xi comes from thermally activated pairs. 
They are given by 

dm 



Xo 

A v (5+ if dwj 
9pe& J?,- Ewi{l + vn) 85 ' { ] 

where we have assumed p+ = pP for simplicity. Figures [| and | show the susceptibilities xo 
and Xi as functions of T and 5 = 3u — 1 (with < 5 < 1/5 corresponding to 1/3 < v < 2/5), 
for three different counting schemes for statistics matrix. Indeed, for 1/3 < v < 2/5, Xi ocT 
and tends to zero as T — >• 0, so that x ~^ Xo(T = 0) as expected. When A ~ 100 mK, 
p e ~ 10 11 cm -2 and ~ p$ (Bohr's magnetic moment), the ratio of the coefficients of xo, 
A = p+/3(f) , and that of Xi, A\ = A/9p e 0Q, is Aq/A\ ~ 10 2 . Therefore, in a certain range 
of filling factor greater than z/ = 1/3, the de Haas- van Alphin term dominates, which is 
positive at T = 0. On the other hand, from Figure we see that Xi is always negative in the 
range of fillings at hand. As T increases, the magnitude of Xi increases, gradually becomes 
comparable to xo an d eventually dominates. Then the total susceptibility will change sign. 
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Thus, we observe that for a given filling factor, there can be a paramagnetism-diamagnetism 
transition, with the critical temperature generically T c ~ A. The new incompressible state 
with z/ = 2/5atT = 0is totally diamagnetic, because with v — > 2/5 one has w^O, which 
causes xi divergent. This might be viewed as a generalization of Landau diamagnetism to 
quasielectrons. 

VIII. CONCLUSIONS AND DISCUSSIONS 

Quasielectrons and quasiholes in the FQH liquids obey fractional (including nontrivial 
mutual) exclusion statistics. Their statistics matrix near the magic electron filling 1/m 
(with m odd) can be determined from various state-counting schemes. These schemes in- 
volve different assumptions on statistical correlations between quasiparticles, resulting in 
somewhat different mutual statistics. A common feature of the schemes is that both quasi- 
holes and quasielectrons in the incompressible FQH liquid background behave like a charge 
in an effective magnetic field. The common assumption is that quasielectrons are in the low- 
est available Landau level (with respect to the effective magnetic field). Then the thermal 
activation of FQH quasiparticles at low temperatures is discussed in the dilute generalized 
ideal gas approximation. If these quasiparticles dominate the low-lying excitation spectrum, 
their contributions dominate the low-temperature thermodynamics. Otherwise, contribu- 
tions from other low-lying elementary excitations (such as skyrmions, if they exist and are 
important) have to be added to the quasiparticle contributions we obtain here. At higher 
temperatures, our assumptions obviously break down. We hope the situation could be im- 
proved in the future by incorporating corrections. Right now we just leave the question of 
when the corrections should be important open to experiments. 

We have used the approximation of generalized ideal gas in our treatment of statistical 
thermodynamics. It is good only when the thermally activated quasiparticle pairs are not 
too dense. On the other hand, to look for the effects of mutual statistics, the densities of 
quasiparticle pairs should not be too low. We hope there is some intermediate range for 
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quasiparticle densities in which the two conflicting requirements could be reconciled to some 
extent. Whether this is true, only experiments can tell. 

It is shown that the existence of hierarchical states in the FQH effect can be viewed 
as a manifestation of the exclusonic nature of the relevant quasiparticles. The associated 
FQH phase transition is shown to occur only at zero temperature. Thermal activation 
of quasiparticle pairs and thermodynamic observables are numerically studied with three 
possible statistical matrices. At zero temperature, they are all equivalent to each other, 
but differences show up at finite temperature. In particular, we have demonstrated that 
for a fixed filling factor between 1/3 and 2/5, with increasing temperature, the system may 
possibly exhibit a transition from paramagnetism to diamagnetism. However, we should be 
cautious about this possibility: We have assumed that /z± are of the same order of magnitude 
as the Bohr magneton; also the approximation of generalized ideal gas might break down at 
the would-be transition temperature. 

It is desirable that these theoretical predictions would be put to experimental tests, if the 
tremendous difficulties in measuring thermodynamic quantities of a thin layer of electron 
gas could be overcome someday. 
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FIGURES 

FIG. 1. The two dimensional plot of the occupation number of the quasielectron n+(T,5). 
The filling factor is near 1/m = 1/3. (a) for the bosonic vortex scheme, (b) for the composite 
fermion scheme, (c) for the correlated vortex or projected composite fermion scheme. In all the 
figures, 5 = 0.15, 0.1, 0.05, and respectively, from above. 

FIG. 2. The two dimensional plot of the occupation number of the quasihole ra_(T, 5). The 
filling factor is near 1/m = 1/3. (a) for the bosonic vortex scheme, (b) for the composite fermion 
scheme, (c) for the correlated vortex or projected composite fermion scheme. In all the figures, 
5 = 0, 0.05, 0.1, and 0.15 respectively, from above. 

FIG. 3. The thermodynamic potential = pV near the filling factor 1/3. (a) for the bosonic 
vortex scheme, (b) for the composite fermion scheme, (c) for the correlated vortex or projected 
composite fermion scheme. In all the figures, 5 = 0.15, 0.1, 0.05, and respectively, from above. 

FIG. 4. The temperature dependence of the specific heat C V (T,5) near the filling factor 1/3. 
The different curves are different <5's. (From above, 5 = 0, 0.05, 0.1, 0.15.) (a) for the bosonic 
vortex scheme, (b) for the composite fermion scheme, (c) for the correlated vortex or projected 
composite fermion scheme. In all the figures, 5 = 0, 0.05, 0.1, and 0.15 respectively, from above. 

FIG. 5. The two dimensional plot of xo(T, 5) (de Haas-van Alphin term of the susceptibility ) 
near the filling factor 1/3. (a) for the bosonic vortex scheme, (b) for the composite fermion scheme, 
(c) for the correlated vortex or projected composite fermion scheme. 

FIG. 6. The two dimensional plot of\;i(T, 5) ( a pair excitation term of the susceptibility ) 
near the filling factor 1/3. (a) for the bosonic vortex scheme, (b) for the composite fermion scheme, 
(c) for the correlated vortex or projected composite fermion scheme. 
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